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INTRODUCTION 


The  Initial  simulation  research  work  done  in  the  MIRADCOM  Infrared 
Simulation  System  (IRSS)  was  heavily  partitioned  between  analog  and  dig- 
tial  computers  that  were  driving  real  time  hardware  [1].  Several  dis- 
advantages to  this  method  of  driving  the  IRSS  hardware  were  noted. 

These  disadvantages  centered  primarily  on  the  complexity  and  maintenance 
of  such  a distributed  system  and  on  the  resulting  control  difficulties. 
It  also  became  apparent  that  complexity,  maintenance,  and  these  control 
difficulties  could  be  significantly  reduced  by  placing  all  computer 
functions  external  to  the  IRSS  in  one  existing  high  speed  digital  com- 
puter that  communicated  with  the  IRSS  facility. 

Only  one  problem  remained  before  this  idea  could  be  Implemented. 

The  airframe  solution  frequencies  were  too  high  for  conventional 
real-time  digital  simulation  methods.  This  report  documents  the  model- 
ing method  that  was  discovered  to  be  suitable  for  solving  this  problem. 

2.  AIRFRAME  EQUATIONS 

The  airframe  equations  originally  developed  for  analog  and  hybrid 
simulation  are  suitable  for  digital  simulation  with  minor  modification. 
These  equations  can  be  reduced  to 

ii'  - + v'r"  - w'q'  - sin  (©L  + 6p)  (1) 

v'  - v'  - £6  r"  + £21  6wi  sin(J>  (2) 

- S5  m'  * q'  ♦ $18  cos  (0L  + 0p)  - £21  cos$  (3) 
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r ' 

■ {1  + « 2 1'  + ?20  6WJ  cos* 

(4) 

| v + C2  r + £2q  6wl  sin({> 

(5) 

:if:;trroiir  "issue  coordi"'“<!  w..  shom  ln  Fi . 1 

"1“°"S  “Sed  ld  (1)  through  (5,  are 

Utility 

t 9 

c pSu"  I * -l 

1 2ly  Cmtotal  + (XCG  - XR)  Cntl 

,1 

(6)  M 

r _ PSu'd  . 

2 21  Cmq 
y w 

1 

(7)  I 

• * : 

y* 

f ^5  " ~ ^ Cn 

3 2m  total 

1 

(8) 

!j 

1 m 

C6  = u' 

r j- 

(9) 

)■  j 

C18  3 8 = 32.174  ft/sec2 

» 

] 

(10) 

j «20  ’ 3'81972  + G2o)  + 6ll(gl  + g2<l)J  ^ 

‘S 

(11) 

^21  3.81972  8c(g  + g «)  T 

x ^ m 

(12) 

8,  - 
1 2 

n 

(13) 

4 

J 

I 


b2  - 0J  (1/Iy) 


3 4 * ~ 6j  ( 1 /in) 


b5  - e4  „• 


b6  - b2  u- 


8e  - 82  «cc  ' V 


B10  ■ e6  d 


611  * 68  u 


«14  ' Th  (1/">  + B5  Ca 


* * 

Cn,,.„f  , - Cm  + Cm 
total  pe 


* * 

Cnt-^.,i  = Cn  + Cn 
total  pe 


Cm  = F 


1 + F2“  + (F3  + F^lfiwilT 


Cn  - f,  + f0u  + (f..  + f.a)  6 .|t 
12  J 4 ' wi1  m 


CD  +“a 
o 


ACa  - (Aj  + A2a)|6wi| 


where  is  a wing  incidence  control  toning  factor  that  must  be  deter- 

mined by  tuning  to  match  test  flight  maneuverability  in  conjunction  with 
TCp  which  is  an  indirect  measure  of  CP  location  that  must  also  be  tuned 


for  matching  test  flight  trim  conditions. 


l£P  is  implemented  in 


^XCG  ” XR^  " (XCG  XR* TABLE  + TCP 


Note  that  Equations  (1)  through  (5)  are  coupled,  nonlinear,  and  time- 
varying.  Additional  equations  for  use  in  the  airframe  model  are 


♦ ■ IP  dt 


= Jq'dt 

a/*-^Vdt 

J cos  0p 


where  p“  is  a roll-rate  time  history  from  a nominal  test  flight. 
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3. 


MODEL  DEVELOPMENT 


Consider  coupled  Equations  (3)  and  (4)  together. 


q'  = £ L w'  + ^2  <l'  + ^20  6wi  cos(*  (32> 

w"  = £5  w"  + C6  q"  + ^18  cos (0^  + 0p)  - S21  6wl  cos<{)  (33) 


i 


Take  the  Laplace  transforms  of  Equations  (32)  and  (33),  and  assume  con- 
stant coefficients  to  obtain. 


q'(s)  “ lq'(o)  + ^ w'(s)  + C2Q  L (6yi  cos<f>)]/(s  - (34) 


and 


w'(s)  = (w'(o)  + f,6  q ”” (s)  + L[518  cos (0^  + 0p) 


- £ 


21  wi 


cos4i] }/ (s  - £,.) 


(35) 


after  rearrangement.  The  coupling  is  shown  in  Figure  2.  Take  the  7- 
transform  of  Equation  (34)  as 


+ ('2oZ[l  (6wt cos^  (r~q)  J 


06) 


Hu*  lirst  2-transform  on  the  ri8ht-hand 
and  can  be  obtained  from  tables,  but  the 
mate  convolution  method  Is  appropriate. 


side  of  Equation  (36)  Is  exact 
last  two  are  not  and  an  approxi- 
The  exact  transform  is 


( 17) 


whe  re 


the  definition  of  the  Z- transform 


|2|  is  taken  as 


Z'OKt)])  - £ f<„T>«"  . f(«) 

n-0  ( 

Note  that  T is  the  simulation  time  step. 

Tunable  trapezoidal  convolution  [21  In  defined  by  Olck.on  «, 
Zf(Os<»)l  * Tf(z>„(z)  - T(,]f  (o)g(z)  e (i  . u)g(0)f(,)| 


7“  b‘'  aWar'  th“'  D1CkSOn  f unable  trapezoidal 

convolution*  Equation  u„,  „ ..EuJer  tun>1)le  convolutlon  „ 

...re  recent  reference  |„.  „e  „1,„  „nen  the  tern,  "tr.pezolde,  

"7 r:"  'TTC,''  th“'  131  “ — * 

not  used  In  this  work. 
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If  Equation  (39)  is  applied  to  the  middle  RHS  term  of  Equation  (3ft), 


- T(l  - n) i. 


(»A) 

V ' t -0 


(*) 


(40) 


Nott'  that 


(41) 


from  the  transform  table  in  the  Appendix.  Equation  (40)  becomes 


zh"’  (’-^)] ' r:!! 


.T?jy,l^pT  _ t(i  - n)w'(i)  . 

*e  1 - re  2 


;'»T 


(42) 


If  Equation  (39)  la  applied  to  the  last  RHS  term  of  Equation  (3ft), 


i 


Z[L<S»t  CO.,)  (.4^]^.,  co.,)j  ^ —KT2T ) 


T( l - n)  i 


OOBlf) 


(43) 


zf 


L(6wi  cos4»)(i-^-)] 


tZ<6w1  cos<t>)  Tn(6„4  cos(}>); 


wi  t-0 


2 fj  V 

l - ze 


1 - ze 


S,T 


T(1  - n)  r/u  (6^  cos4>) 


(44  ) 


Substitute  Equations  (37),  (42),  and  (44)  into  Equation  (36). 


*'M  * h 

1 - ze 


Tw"(z)  TQw^o) 

V " ^2T 

1 - ze  1 1 - ze  2 


- t(i  - n)w 


-(.)] 


+ e 


20 


TZ<6wi  cos<J>)  Tn(6ttl  cost})) 


1 - 


ze 


wi__^t=0  , ry 

- T(1  - n)2(6wl  cos0) 

1 - ze 


(45) 


> £ T 

_ q"(o)  + ^1T[w"(z)  - nw'(o)  - (1  - n) (1  - ze’2  ) w'(z)| 
q (z)  - — - — ■ — 


1 - ze 


+ ?’20T[^(6wi  C°S><,)  ~ n(6wi  coa<t>)t-0  ~ (1  ~ 0)  (l~ze  2 )Z<fiwt«-^)] 


f t 

1 - ze  2 


(46) 
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From  Equation  (38),  it  is  seen  that 


w'(z) 


and 


%x> 


l 

n-0 


w"(nT) zn 


(47) 


Z(6wi  COS(*,)  “ I 


n-() 


l^i  tos(ionT)]zn 


(48) 


and 


q'(z)  - l q ”*(nT)  zn  . 
n-0 


(49) 


Substitute  Equations  (47),  (48),  and  (49)  Into  Equation  (46);* 


(C  t\  00 

1 " ZC  2 /n-0  ^(nT)  Z"  “ ^(0)  + ^ f t w>(''T)2«  - nv'(o) 

' u Ln-o 


- (l  - n) 


(‘ ' “It  ”'<nT)'n] 


+ ^20TJ_X„  (nT)  cos(o)nT)  ) zn 


n-0 


n(6wl  cos^)t>0 
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(i  - n) 


(■ 


S,t\  " I 

ze  I I (nT)  cos(umT)]  z | 

/ n-0  ! 


(50) 


cx> 


I 

n-0 


q'(nT)zn 


n=0 


q^(nT) 


n+1 

z 


- q'(o)  + CxTtl  - (1  - n)]  I w'*(nT)z11  - C1Tnw'(o) 

n=0  1 


£ T 00 

+ CjTd  - n)  e 2 I w'(nT)zn+1 
n =0 


+ S20T[1  - (1  - n)] 


00 


I 

n-0 


[6  (nT)  cos  (umT)]zn 

W1 


■ ^20Tn(<Swi  cos<*>)t=0 


?2t  °° 

+ £20T(1  " n)e  E 1 6 <(nT)  cos(u>nT)]zn+1  (51) 

n=0  W1 

Now  the  Z-transform  shifting  theorem  consistent  with  the  definition  of 
Equation  (38)  is  - • 
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(52) 


in  “ m- 1 

* f(z)  - 7 f(nT  - mT )z"  - J f(nT  - mT)  z'1 

n-0  n-0 


where 


m > 0 


(53) 


and  m is  an  integer.  The  shifting  theorem.  Equation  (52),  can  also  be 
written  as 


OO 


I 

n-0 


f (nT) 


n+m 

z 


l f(nT  - mT)  z" 
n-ra 


(54) 


which  is  directly  useful  in  many  cases.  Use  the  shifting  theorem. 
Equation  (54),  in  Equation  (51). 


OO 


l 

n-0 


q'(nT) 


n 


z 


C,T  <*> 

e l q>(nT  - T)  zn 
n*l 


“ q'(o)  - ^'(o)  - ^20Tn(6wJ  cos^)t-0 

00  £ T oo 

+ CjTn  l w'(nT)  zn  + C T(l-n)  e’2  Y „'<nT  - T)  z" 
n-0  n-l 

OO 

+ C20Tn  l [6wl  (nT)  cos  (u>nT)  ] z" 
n-0 
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> ‘2 


(55) 


+ 40T(1  - n)  e l (6  (nT  - T)  c 
n=l  wl 


os  [u>(nT  - T ))}2n 


Compiete  all  Z-transforms  in  Equation  (55) 


[nL  q ^ " G 2 J0  q'(nT  - T)  ^ q'(-T)] 


^C1w'(o 


-|q  (°)  - m 5^(0)  + c20  (6wl  cos^) 


t-0  J 


+ qTn  n=0  W^<nT)  Z"  + V(1  " n)  e^l  w'(nT  . T)  zn 

n=0 


1 ^ W + ^20^  E cos  CoonT)  J zn 


+ ^20T(1  ~ H)  e 2 l {6vi(nT  - T)  cos  [u)(nT  - T)}zn 


C T 

- £20T(1  - n)  e 2 (6wl(-T)  cos(-u)T)]  j 


The  ability  of  Equation  (56)  to  handle 
be  shown  by  letting 


a nonzero  initial  condition 


n =>  0 


and  noting  that 

q (o)  * q'(o) 
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(58) 


results  when  coefficients  of  z"  are  equated. 


Now  let 


n « 1 


and  equate  coefficients  of  zn  in  Equation  (56)  to  obtain 


^2^  £ T 

q'(T)  * e ‘ q'(o)  + qT(l  - n)  e’2  w'(o) 


+ ^20T^  ~ I <5.  (°)  cos(o)J 


+ CjTnw'd)  + S20Tn[6wl(T)  cos(a)T)  ] . 


If  one  lets 


n = 2 


and  equates  coefficients  of  zn  in  Equation  (56),  then 


q'(2T)  » e ‘ q'(T)  + ^1(1  - n)  e’2  W'(T) 


v^2  * 

+ C2QT(1  - n)e  I 6wi<T>  cos  (WT))  + SjTnw'(2T) 


+ C20TnI6wi(2T)  cos  (2u,T>)  • 


Comparison  of  Equations  (60)  and  (62)  shows  that  for 


n > 1 


(63) 


the  coupled  general  recurrence  equation  is 


qn  = e qVl  + VnV  + V(1  "n)e  Vl 


+ C20Tn(6wi  C08*>n  + *>20T(1  “ n)e  (<5wl  cos^>n-l  * «>4> 


If  an  identical  procedure  is  used  on  Equation  (35),  one  obtains 
the  other  coupled  general  recurrence  equation. 


Wn  = 6 ^ Vl  + Vnqn  + V(1  " Vl  • 

+ Tn(S18  cos(0L  + 0p)  - ^216w1  cos01n 

c5t 

+ T(1  - n)e  [C18  cos(0L  + 0p)  - C216w1  cos*]^  (65) 


for 


n > 1 


(66) 


Also,  if 


n - 0 


(67) 
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the  result  is 


w (o)  • v'(o) 


(68 


for  handling  nonzero  initial  conditions. 


Equations  (64)  and  (65)  are  algebraic  equations  that  can  be  solved 


simultaneously  for  q'  and  w^.  The  results  are 


r v 

- + c 

%l~T- 


!V  (1  - n)ne 


2 2 


W n 


’n-l 


v[ne5  + a - n)e^T]  A,C1ST2  n2  \ 

Vi  + v - wv)  cos(0l 


1 - n‘ 


+ Vn 


C T " 

ci^18t2  (l~n)ne  5 

+ 5 — ? cos(0  + 0 ) 

1 - ZjZfT2  n2  J L p "-1 


, V'  ^ 


fTn^20  - ^1^21Tl1^1 

n - ^ J 


(6  . cos*!1) 

wi  n 


+ [T(1  " n)(j 


C2T  , _ _ ”5 


hoe  ~ ^i^21Tne 


SCT 


2 2 


H 


<6  cos*)  , (69) 


r 


■ 


x , .,fi8T(i-n)/'5'r 

c4  

1_,ci  VV 


L ‘-^frV 


(81) 


c5 


(82) 


the 


TC6  ■ - «<T 

L ^lVV 

aPproximati0ns 


(83) 


c°s(0L  + e ) a 


P'n  - cos(0L  + 0p) 


n-1 


cosC9l  + ep)  . 

P n-l  ~ cos(0  + e ) 

L P^n -2 


are 

(69) 


made  and  Equation. 
and  (70),  then 

^ * b 


^72^  thr°ugh 


are  substituted 


(84) 

(85) 

lnto  Equations 


n ~ + bnW-  + 

0 n-l  + a2c°s (0  + q \ 

P n “I  + alCOs(0r  + e ) 

L P^n-2 


°(6wi  cos0)  + e 

• n + e3(6  cos0) 

and  1 n-l 


(86) 


V ~ 


n e2wn-l  + eiq^  + 


eo  c°s(e,  + e v 

L r'n-l  + T 


+ W<5 


C5V°W  COS0)n  + 


c4  c°s(0  + 0 >. 

1 Vn -2 


n c6(<5wI  cos0) 


n-l 
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(87) 


Equations  (86)  and  (87)  can  each  be  evaluated  with  one  line  of  digital 
program  code. 

Consider  coupled  Equations  (2)  and  (5)  together. 


r = ~€jV  + ^2r  + ^20  ^wi  s*n<^ 


v'  = - Z6r'  + C21  <Sw.  sincf. 


Take  the  Laplace  transforms  of  equations  (88)  and  (89)  and  assume 
constant  coefficients  to  obtain, 

r'(s)  = [r'(o)  - Cjv'(s)  + C2oL^wi  sin<(>)  ] / (s-C2) 


v'(s)  = [v'(o)  - 56r'(s)  + ^21L(6w1  sin<t,)]  / (s"^5)  (91) 

after  rearrangement.  The  coupling  is  shown  in  Figure  3.  By  direct 
analogy  with  the  method  used  on  the  q'(s)  and  w'(s)  equations,  one  can 
find  that 

C T £ t 

r n ~ e Vl  - 2 K-l  + koT^Swi  sin*>n 


+ C20T(l-n)e  (6wl  sin<)))n_1 


. *5T  - 

V'  a e V 


n = e Vl  “ 56Tnrn  " ?6T(1-ri)e  rn-l  + 521Tn<6wl  8in<»n 


+ C21T(l-n)e  («  sin*)^ 
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Equations  (92)  and  (9J)  are  algebraic  equations  that  can  be  solved 
simultaneously  for  t ^ and  v J.  The  results  are 


( T l 

• + '.,f.(lT‘(l-u)n‘’ 

r"  [ 1 - C.^tV 


C,T|noY  ♦ u-n)^  v. 
1 - 


"Tn^20_^lf’21T,,)" 

- '-‘iVV  - 


W«1  sl"*\, 


T(l-n)(f.20e  - i.jt.^jTne 

1 - W2"2 


r C5t  2 C2T1 

e + £ £ T"(l-r()q£  *■ 

"1  1-WV  JVl 


(6wi  8ln^)n-l 


C2t  £.  T 

C6T|ne  +(l-n)e  ] 


1 - CiVV 


Tti(t,21  - i,6t,20Tn> 


+ l 5 I (<5  . slntf) 

I 1-^jV  I wl 


T(i-n)q21o  -£6C20Tne  ‘ ) 

1 - W2n2 


(6wl  sin*Vl 


n > 1 


(100) 


Substitute  Equations  (72),  (73),  (76),  (77),  (78),  (79),  (82),  and  (83) 
into  Equations  (98)  and  (99).  The  results  are 


rn  = blrn-l  ~ boVn-l  + ao(6wi  8in(*>n  + e3(<Swi  sin(<>)n-l 


(101) 


and 


vn  = e2Vn-l  - elrn-l  ~ Tc5<6wi  sin^n  “ Tc6(6wl  sln0>n-l  • <102) 


Equations  (101)  and  (102)  can  each  be  evaluated  with  one  line  of  digital 
program  code. 


If  Equation  (30)  is  written  with  an  initial  condition. 


0p  - f q'dt  + eSE  . 


(103) 


Take  the  Laplace  transform  of  Equation  (103), 


q"(s)  + 

ep(»)  ■ i—56 


(104) 


Take  the  Z-transform  of  Equation  (104) 


eP(*)*Z[q'(»)(i)]  + eSEZ4> 


(105) 


Use  Equation  (39)  and  the  transform  table  in  the  Appendix  on  Equation 
(105). 


6p(z)  = Tq  (z)  (j“)  - TlnqtoMj^-) 


+ (l-n)L-1(|‘)t=oq"(z)]  +I!f 


(106) 
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(107) 


9p(z)  » 


Tq'(z)  - Tnq'(o)  - T( 1-n) ( l-z)q(z)  + 0g 
1 - z 


(l-z)0p(z)  a Triq'(z)  + Kl-nWz)  - Tnq'(o)  + 0J 


Applicable  Z-transforms  definitions  are 


9p(z)  “ l ep(nT)zn 

n=0 


(108) 


(109) 


q"(z)  “ I q (nT)zn 
n“0 


Substitute  Equations  (109)  and  (110)  into  Equation  (108) 


(l-z)E  Op(nT)zn  * Tn£  q'(nT)zn  + T(l-n)z£  q'(nT)zn 


- Tnq(o)  + 0, 


(110) 


(111) 


I 0p(nT)zn  - l 0 p(nT)zn+1  « Tnl  q'(nT)zn  + T(l-n )£  q'(nT)zr 
n“0  nm0  n“0  n-0 


- Tnq(o)  + 0_ 


Use  the  shifting  theorem.  Equation  (54),  in  Equation  (112). 


(112) 


00  OO  oil 

1 ep(nT)zn  - l 0„(nT-T)zn  - Tn£  q'(nT)zn 

n-0  n-1  n=0 


+ T(l-n)I  q'(nT-T)2n  - Tnq^o) 
n=  1 


+ 6SE  <113> 


Complete  all  Z-transforms  in  Equation  (113). 

00  OO  OO 

1 0p(nT)zn  - l e (nT-T)z11  + B„(-T)  = ml  q'(nT)zn 

n-0  n=0  1 n-0 


+ T(l-n)I  q'(nT-T)zn 
n*0 


- T(l-n)qt-T)  - Tnq'Co) 

+ eSE  (H4) 

The  ability  of  Equation  (114)  to  handle  a nonzero  initial  condition  can 
be  shown  by  letting 


n-0  (115) 

and  noting  that 

0p(o)  = 0gE  , (superelevation  angle)  (lib) 

results  when  coefficients  of  zn  are  equated.  Now  let 

n - 1 (117) 

and  equate  coefficients  of  z1'  in  Equation  (114)  to  obtain, 

6p(T)  - 0p(o)  + Tnq'(T)  + T(l-n)qto)  . (118) 
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Now  let 


n - 2 (119) 

and  equate  coefficients  of  zn  in  Equation  (114)  to  obtain, 

Bp(2T)  = 0p(T)  + Tnq'UT)  + T(l-n)q'(T)  . (120) 

Comparison  of  Equations  (118)  and  (120)  shows  that  for 

n > 1 (121) 

the  general  recurrence  equation  is 

6p  ~ 0p  + Tnq'n  + T(l-n)q'_ x . (122) 

n n-1 


Equation  (31)  can  be  written  with  an  initial  condition  as 


/ x9  dt  + ^lead 
o 


(123) 


where 


T9  - r'/cos9p  . 


(124) 


Take  the  Laplace  transform  of  Equation  (123). 
19(3)  + ^lead 

Vs)  — — 


(125) 


If  Equation  (125)  is  treated  the  same  as  Equation  (104),  then 


*P  = *p  + TnT^  + T(l-n)T9n.j 
n n-1 


(126) 


and 


V°>  * ^LEAD’  ^lead  angle) 


(127) 
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Make  the  definitions. 


Tcl  - Tn  (128) 

and 

Tc2  = T(1-ri>  ’ (129) 

Substitution  of  Equations  (128)  and  (129)  into  Equations  (122)  and  (126) 
yields 


Tcl9n+  Tc29n-1 


and 


(130) 


% ’ ‘‘'p  . + Tcl  T9n  + tc2  T9n-1 

n n-1 

for 

n > 1 . 


(131) 


(132) 


Equations  (130)  and  (131)  can  each  be  evaluated  by  one  line  of  digital 
program  code. 


Now  the  longitudinal  acceleration  of  the  missile  can  be  computed 
algebraically  with  the  aid  of  Equations  (1),  (21),  (86),  (87),  (101), 
(102),  and  the  approximation 

sin(0L  + 0p)n  « sin(0L  + Qp)^  (133) 


as 


“n  = C14  + Vn  " V>n  “ C18  sin(0L  + 6P>n-l  * 
Tunable  integration  [2]  of  Equation  (134)  yields 
un  = un-l  + + 


(134) 


(133) 


30 


where 


u'(o) 


UBORECLEAR  ' 


(136) 


Note  that  the  tunable  Integrator  is  in  the  same  form  as  that  obtained 
by  tunable  trapezoidal  convolution  (2]  of  integration  of  uncoupled 
variables.  They  are  equivalent.  Substitute  Equations  (128)  and  (129) 
into  equation  (135). 

U - u . + T , u + t ,,u  , ^117^ 

n n-1  cl  n c2  n-1  viJ/7 

Equations  (134)  and  (137)  can  each  be  evaluated  with  one  line  of  digital 
program  code. 

Finally,  Equation  (29)  can  be  written  with  an  initial  condition  as 
t 

^ / P dt  + ^BORECLEAR  (138) 

o 

and  treated  the  same  as  equation  (103).  Usually, 


^BORECLEAR  “ 0 


(139) 


The  resulting  recurrence  Equation  for 
n > 1 

is  then 


(140) 


♦n  * *n-l  + TclPn  + Tc2Pn-l  d*l) 

where 

" ^BORECLEAR  * d*2) 

The  p roll  rate  history  to  be  integrated  must  be  representative  of  test 
flight  conditions  and  be  available  over  the  missile  flight  time  of 
interest. 
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4. 


DISCUSSION 


Although  the  airframe  model  has  been  derived  with  the  assumption  of 
constant  coefficients,  it  is  known  that  the  coefficients  are  time  varying 
and  that  some  have  discrete  jumps  corresponding  to  wing  erection  and 
motor  events.  In  order  to  minimize  the  error  associated  with  such  an 
assumption,  the  coefficients  are  updated  only  one-fourth  as  often  as  the 
airframe  variables.  This  method  allows  the  coefficients  to  appear  con- 
stant for  three  out  of  every  four  passes  through  the  airframe  equations, 
but  yet  allows  time  variation  and  thus  discrete  jumps  in  these  coefficients 
on  the  fourth  pass.  This  approach  is  desirable  not  only  because  it  allows 
stable  solutions  when  the  tunable  trapezoidal  convolution  [2]  phase 
parameter,  n,  is  tuned,  but  also  because  it  relieves  the  real  time 
computing  load  by  not  requiring  coefficient  recomputation  in  each  pass 
through  the  airframe.  Passes  are  identified  with  frame  numbers  that 
cycle  through  one  to  four  as  shown  in  Figure  4.  The  coefficient  comput- 
ing load  is  then  distributed  among  the  four  frames.  This  method  also 
allows  the  distribution  of  other  simulation  computations  in  the  four 
frames  as  long  as  real-time  computing  limits  are  not  exceeded. 

This  airframe  model  has  been  used  in  three  applications  to  date  as 
illustrated  in  Figure  5.  The  first  application  was  a real  time  hardware- 
in-the-loop  IRSS  simulation.  This  simulation  was  implemented  as  shown 
in  Figures  5a  and  6.  The  purpose  of  the  airframe  model  is  to  drive 
IRSS  commands  that  cause  real  time  motion  of  the  hardware-in-the-loop 
(guidance  section).  The  IRSS  then  supplies  measured  hardware  motion  as 
well  as  guidance  section  output,  5wi  .,  in  order  to  close  the  simulation 
loop.  It  should  be  noted  that  the  guidance  section  is  responding  to  a 
real  time  projected  IR  environment  during  this  time  and  is  therefore 
driving  the  simulation  when  in  closed  loop.  Note  also  that  open  loop 
capability  is  present  for  checkout  and  other  studies. 
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The  second  application  is  to  emulate  the  IRSS  system  as  shown  in 
Figure  5b.  The  emulation  is  critical  to  use  of  the  IRSS  simulation  since 
it  allows  checkout  and  some  validation  to  be  accomplished  without  sub- 
jecting the  actual  guidance  section  to  heavy  wear  during  checkout.  Also, 
a majority  of  the  IRSS  checkout  can  be  emulated  even  when  the  IRSS  facil- 
ity [1)  has  frequent  maintenance  or  other  missile  users.  Then  all 
lessons  learned  on  the  emulator  are  transferred  to  the  IRSS  program  with 
its  identical  airframe. 

The  third  application  is  to  simulate  a rolling  airframe  missile  in 
a pure  digital  simulation  as  shown  in  Figure  5c.  An  effective  digital 
seeker  model  is  required,  but  extra  computation  spent  here  is  offset  by 
reduced  computation  in  the  airframe  model  presented.  The  result  is  a 
reasonably  economical  simulation  suitable  for  kinematic  boundary  studies 
and  other  general  applications.  It  should  be  noted  that  the  seeker  model 
used  here  is  also  used  in  the  IRSS  emulator. 


Figure  A.  Airframe  model  calculation  block  diagram. 
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IRSS  Emulator 


EXTRA  WLA TOR  i 1MIT 

U |4wln!  ^wi'uMlT 
THEN 

iw‘n  ’ Awi  LIMIT 


/'nail,  i ■ . i 


<6  . co*4)  , 

wl  n- 1 


"‘n-l  { MEASURED  FROM  HARDWARE 

MEAS  j 

> EXTRAPOLATED  WING  INCIDENCE 

I <tn  - M-Swl,,.,  -ivl„.2)  ♦ «vin.3 


4win-3  - 6wln.2 
4v,1n-2  " 4“Vl 


CLOSED 

LOOP 


ROLL  POSITION  ♦ COMMAND 
'n  *n-l  + Tclpn  + 1c2pn-l 

r 

[ REDUCE  TO  PRINCIPAL  ANCLE 


Figure  6. 


Airframe  model  implementation  in  the  IRSS.  (Continued) 
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SYMBOLS 


A1‘  A2 


AC. 


DO 

i i 

Cm,  Cn 


Cmw,  Cnw 
Cm 

q 

* * 

Cm,  Cn 


Cm  , Cn 
pe  pe 


* * 

Cin  „ - , Cn 
total 


fr  f2’  f 3*  f4 


Fl*  F2’  F3’  F4 


8^  * ^2 


<V  G2 


Curve  fit  functions  for  representing  AC. 

A 

Axial  drag  coefficient 

Incremental  axial  force  coefficient 
Zero  lift  drag  coefficient 

Moment  and  normal  force  coefficients 

Control  wing  moment  and  normal  force  coefficients 

Damping  coefficient 

Modified  secant  slope  moment  and  normal  force 
coefficients  without  plume  effects 

Exhaust  plume  effect  moment  and  normal  force 
coefficient  increments  to  c£  and  Cn 

Total  modified  secant  slope  moment  and  normal 
force  coefficients 

Missile  reference  diameter 

Base  of  natural  logarithms 

Curve  fit  functions  for  representing  Cn 

Curve  fit  functions  for  representing  Cm 

Acceleration  due  to  gravity 

Curve  fit  functions  for  representing  Cnw^ 

Curve  fit  functions  for  representing  Cmw^ 

Moment  of  Inertia  about  y,  or  z axes 

Missile  roll,  pitch,  and  yaw  rates 


M 


IRKIUNO  PiOS  Mi.Aiy 
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s 

t 


U , 

V , 

w 

• 

• y 

• 

U , 

V , 

w 

x,  Y,  Z 


*r 


X',  Y",  Z' 


a 


V 


P 


P 


<t>' 


( ) 
( ) 


6=0 


t=0 


Missile  reference  area 

Time  from  missile  boreclear 

Missile  thrust 

Missile  velocity  components 

Missile  acceleration  components 

Missile  fixed  coordinates 

Distance  between  CG  and  missile  nose 

Distance  between  wind  tunnel  reference  point  and 
missile  nose 

Nonrolling  missile  coordinates 

Missile  angle  of  attack 

Missile  control  wing  incidence 

Euler  angles  from  earth  fixed  to  launch  coordi- 
nate system 

Euler  angles  from  nonrolling  missile  to  launch 
coordinates 

Atmospheric  density 

Roll  angle  between  missile  fixed  and  nonrolling 
missile  coordinate  axes 

Roll  angle  between  missile  fixed  coordinate 
system  and  the  projection  of  the  missile  velocity 
vector  in  the  Y^Z"  plane 

At  zero  wing  incidence 

At  time  from  boreclear  equal  zero 
Coefficients 
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V Tc4 


i 


< >n 

At  present  time 

< >n-l 

At  present  time  minus  T 

< >n-2 

At  present  time  minus, 2T 

< >n-3 

At  present  time  minus  3T 

T 

Simulation  time  step 

( ^ BORECLEAR 

At  time  t « 0 

n 

Tuning  parameter 

n 

Integer  time  step  number 

(n  - 0 @ t ■ 0) 

0SE 

Superelevation  angle 

\eAD 

Lead  angle 

z 

Complex  variable  used  in 

Z-transform  definition 

s 

T~9 

Complex  variable  used  in 
definition 

laplace  transform 

Z 

Denotes  the  operation  of 

taking  a Z-transform 

L 

Denotes  the  operation  of 
transform 

taking  a laplace 

L*1 

Denotes  the  operation  of 
laplace  transform 

taking  an  inverse 

Cj.  s2*  s5>  s6. 

Utility  variables  and  constants 

^18*  ^20*  ^21*  Bl» 

B2»  B4.  B5,  B6,  Bg, 

610*  Bll'  ^14 
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